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The molecular polar disc in NGC 2768* 
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ABSTRACT 

We present CO(l-O) and CO(2-l) maps of the molecular polar disc in the elliptical 
galaxy NGC 2768 obtained at the IRAM Plateau de Bure Interferometer. The maps 
have a resolution of 2"6 x 2"3 and l'/2 x l'/2 for the CO(l-O) and CO(2-l) lines, 
respectively. The CO maps complete the unique picture of the interstellar medium 
(ISM) of NGC 2768; the dust, molecular gas, ionised gas and neutral hydrogen (H i) 
trace the recent acquisition of cold and cool gas over two orders of magnitude in 
radii (and much more in density). In agreement with the other ISM components, 
the CO distribution extends nearly perpendicularly to the photometric major axis 
of the galaxy. Velocity maps of the CO show a rotating polar disc or ring in the 
inner kiloparsec. This cool gas could lead to kinematic substructure formation within 
NGC 2768. However, the stellar velocity field and H/3 absorption linestrength maps 
from the optical integral-field spectrograph SAURON give no indication of a young 
and dynamically cold stellar population coincident with the molecular polar disc. Very 
recent or weak star formation, undetectable in linestrengths, nevertheless remains a 
possibility and could be at the origin of some of the ionised gas observed. Millimetre 
continuum emission was also detected in NGC 2768, now one of only a few low- 
luminosity active galactic nuclei with observed millimetre continuum emission. 

Key words: galaxies: individual: NGC2768 - galaxies: elliptical and lenticular, cD 
- galaxies: ISM - galaxies: evolution - galaxies: kinematics and dynamics - galaxies: 
structure 



1 INTRODUCTION 

Galaxy formation scenarios must explain the rather uniform 
global properties o f early-type galaxies as seen in the funda- 
ment al plane (e.g. iDressler et al.lll987l : iDiorgovski fc Davis] 
1 19871 ) while simultaneously accounting for the diversity 
seen in the internal k inematics and stellar populations (e.g. 
Ide Zeeuw et alJl2002h . Indeed, the maps of the stellar kine- 
matics of 48 elliptical and SO galaxies from the SAURON 
optical integral-field unit (IFU) have shown a wide vari- 
ety of kinematic substructures, from large embedded discs 
to small counter-rotat ing cores (e.g. lEmsellem et al.l 12004 : 
iMcDermid et al.ll2006h . Major mergers have been found to 
reproduce many of the kinematic substructures observed in 
the SA URON maps witho ut invoking post-merger star for- 
mation (jjesseit et al.ll2007h . However, many of the kinemat- 
ically decoupled c omponents (KDCs) in the SAURON maps 
are very young l|Kuntschner et all 20061 : IMcDermid et all 
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2006), suggesting their stars formed in situ, and not as part 
of another galaxy. 

Over the past twenty years, it has been shown that 
E/SOs do have detectable amounts of cool and cold gas. Cold 
dust was first found in the IRAS survey in which around 45% 
of ellipticals were de tected in both the 60 and 100 /J.m bands 
jKnapp at all 1 19891 ). Molecular gas detections have often 
been biased toward dust-rich or infrared-bright galaxies, but 
two recent surveys do not select on these criteria and give 
detection rate s of 28% for E/SOs in the SAUR ON represen- 
tative sample (jCombes. Young fc Bureaull2007 ) and 7 8% for 
SOs in a volume- limited sample I Sage fc Welch! [2006). Par- 
ticularly interesting for substructure formation are the maps 
of cold molecu lar gas provided b y interferometric CO ob- 
servations (e.g. lYoundl2002j , 120051 ). These maps often show 
rotating molecular gas discs on the same spatial scales as 
KDCs. 

While the current evidence indirectly connects molecu- 
lar gas discs with young KDCs, the details of the connection 
are still unknown. Do all molecular discs form stellar com- 
ponents? Do stars form everywhere in the molecular disc, 
or only in the densest regions? Does the regularity of the 
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Table 1. Basic properties of NGC 2768. The left and middle 
columns list the different quantities and their values; the right 
column lists corresponding references. 



Quantity 


Value 


Ref. 


R.A. (J2000.0) 


09 11 37.413 


1 


Dec. (J2000.0) 


+60 02 14.86 


1 


Heliocentric Velocity 


1373 ± 5 km s" 1 


2 


Distance 


21.5 Mpc 


3 


Scale 


1"= 104 pc 


3 


Type (RC3) 


E6 


4 


Type (CAG) 


SO 


5 


Corrected apparent B mag 


10.60 


6 


Corrected absolute B mag 


-21.15 


6 


L B 


2.3 x 10 10 L Q 


7 


Lfir 


4.3 x 10 s L 


7 


Lfir/Lb 


1.9 x 10~ 2 


7 


Lfir/Mh 2 


6.3 L /M Q 


7 



Refere nces: (TllNagar et al.N2005|) ; (2) N ED; C 3) lEmsellem et al, 
20041') ; (4) Ide Vaucouleurs et al.l Jl99lf) ; (5) ISandage fc Bedke 



19941 ); (6) HyperLEDA; (7) D erived quantity using data from 



NED, IRAS jMoshir et a,l.|[l99 0h. and this paper 



morphology of the molecular disc (presumably indicating 
how long since the gas was accreted) correlate with the age 
of the stars? To answer these questions, we must compare 
the molecular gas morphology and kinematics to the stel- 
lar morphology, kinematics and ages in a spatially-resolved 
manner. This comparison requires interferometric CO maps 
and optical (or near-infrared) IFU data for many individual 
galaxies. With this comparison in mind, we have started to 
obtain CO synthesis maps for the 12 early-type SAURON 
galaxies recently detected in CO (see ICombes et al.l 120071 . 
and references therein). 

In this paper, we present the first CO synthesis maps of 
NGC 2768 and compare the molecular gas to the stellar and 
ionised gas properties from S AURON. NGC 2768 is class i- 
fied as an E6 galaxy in RC3 llde Vaucouleurs et al.lll99ll) 
an SO in the Carnegie Atlas of Galaxies (ISandage fc Bedkd 
1 19941 ). and as a fast rotator i n the recently proposed kine- 
matic classification scheme of lEmsellem et alj (|2007t ) . iKiml 
(1989) first discovered both the polar orientation of the in- 
ner dust lane and the rotation of the ionized gas about 
the major axis, suggesting an external accretion origin 
for the dust and gas. NGC 2768 has a low-luminosity 
active galactic nucleus (AGN) with a low-ionis ation nu- 
clear emission-line region ( LINER) spectrum (iHeckmanl 
ll980T l. a compact radio core (|Nagar. Falcke fc Wilsonll2005l ) 
and an X-ray source consistent with being a point source 
l|Komossa. Bohringer fc Huchral Il999l ). General properties 
of NGC 2768 are listed in Tabled] 

This paper is organized as follows. In Section 2, we de- 
scribe the calibration and reduction of the data from the 
Plateau de Bure Interferometer. Section 3 presents flux mea- 
surements for the continuum data and channel maps, total 
intensity maps and velocity fields for the CO line data. In 
Section 4, we discuss the origin and fate of the molecular 
gas and also consider the contributions to the mm contin- 
uum flux from the low-luminosity AGN (LLAGN) and cold 
dust. We draw our conclusions in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed NGC 2768 in December 2005 (C configura- 
tion) and August 2006 (D configuration) at the IRAM 
Plateau de Bure Interferometer (PdBI). The 12 CO(1-0) and 
12 CO(2-l) lines were observed simultaneously using dual- 
band SIS receivers. All six 15m antennae were used for the 
C-configuration observations in December while only five an- 
tennae were available for the August observations. Typical 
system temperatures were 200-300 K for the CO (1-0) line 
in both months and the CO(2-l) in December. In the worse 
weather of August, the CO(2-l) line observations had sys- 
tem temperatures of 300-1000 K. The spectral correlators 
were centred at 114.743 GHz and 229.482 GHz, the transi- 
tion frequencies of CO(1-0) and CO(2-l) corrected for the 
galaxy's heliocentric velocity of 1373 km s . The correla- 
tor configuration used four 160 MHz-wide (417 and 209 km 
s" 1 for CO(1-0) and CO(2-l), respectively) units with a fre- 
quency resolution of 1.25 MHz (3.3 and 1.6 km s _1 ), covering 
a total bandwidth of 580 MHz (1670 and 835 km s" 1 ). The 
correlator was regularly calibrated by a noise source inserted 
in the IF system. 

We obtained visibilities with series of twenty 1 min 
integrations on source, followed by three 45 s phase and 
amplitude calibrations on a nearby quasar, 0836+710. As 
0836+710 was not bright enough, we used 3C84 (25 Decem- 
ber), 0528+134 (25 August) and MWC349 (28 August) as 
flux calibrators. Comparing the flux measured for 0836+710 
on 28 August to that measured on 25 December suggests an 
uncertainty in our flux calibration of « 20%. This is unfortu- 
nately not a very good estimate as the true flux of 0836+710 
may have varied over the 8 month interval. 

The data were reduced with the Grenoble Image and 
Line Data Analysis Sys tem (GILDAS) software p ackages 
CLIC and MAPPING (jGuilloteau fc Lucasj [2000 ). Using 
CLIC, we first calibrated the data using the standard 
pipeline. Data obtained during periods of bad weather were 
flagged and then ignored. After calibration, we used MAP- 
PING to create two data cubes for both lines with velocity 
planes separated by 40 km s _1 and 60 km s . The primary 
beam size is 44" for the CO (1-0) observations and 22" for the 
CO (2-1) observations. We choose the spatial dimensions of 
the datacube to be twice the diameter of the primary beam, 
88" x 88" and 44" x 44" for the CO(1-0) and C(2-l) lines, 
respectively. The synthesized beam sizes are 276 x 273 for 
CO(1-0) and 1"2 x L"2 for CO(2-l). We thus choose spatial 
pixels of .7" x .7" and .35" x .35" for the CO(1-0) and CO(2- 
1) lines, respectively. As we detected continuum emission, we 
subtracted the continuum using the task contsub in the Mul- 
tichannel Image Reconstruction, I mage Analysis and Dis- 
play (MIRIAD) software package (jSault. Teuben fc Wright! 
Il995h . The dirty beam has small side lobes that necessi- 
tated cleaning the da tacube. The cle aning was done using 
the Hogbom method (|Hogboml 1 1974h : we stopped cleaning 
in each velocity plane after the brightest residual pixel had 
a value lower than the rms of the uncleaned datacube. 

To constrain the continuum emission, we selected chan- 
nels at least 40 km s" 1 away from the highest and lowest 
velocity channels with any line emission in the cleaned data 
cubes. The very edges of the 580 MHz bandwidth were also 
avoided. The resulting continuum windows were 356 and 
200 MHz wide at 115 and 230 GHz, respectively. We used 
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Figure 1. 3 mm continuum contours (white) over a Hubble Space 
Telescope (HST) Advanced Camera for Surveys (ACS) image 
(F8f4W filter, black contours) of NGC 2768. Contours are at 2, 
4, and 6cr, where a = 0.28 mjy beam -1 . The synthesized beam 
is plotted in the lower left corner. 



the MAPPING tasks with the same parameters as for the 
line data to create clean images of the continuum at both 
frequencies. We fit point-source models in the uv plane to 
find the flux of the 3 mm continuum and an upper limit for 
the 1 mm continuum. The continuum emission is spatially 
unresolved at both frequencies. 

3 RESULTS 

We detected extended emission in both the CO (1-0) and 
CO(2-l) lines and point-source continuum emission at 3 mm. 

3.1 Continuum emission 

The continuum maps at 1 and 3 mm both show peaks at 
the centre of NGC 2768. The 3 mm continuum map has a 
peak of over 6a consistent with being a point source (see 
Fig. [TJ . Fitting a point-source model in the uw-plane yields 
a flux density of 1.94 mjy at a right ascension of 09 11 
37.38 and declinatio n of +60 02 14.8. T his position is off the 
VLA coordinates of lNagar et all (|2005l ) by -0.5" in right as- 
cension and -0.1" in declination. Given the 2" 6 x 2" 3 beam, 
this source is however still identifiable with the radio source, 
despite the large offset in right ascension. The 1 mm contin- 
uum map shows a weak peak at the same location. Fitting a 
point-source model in the uu-plane to the same coordinates 
as found for the 3 mm continuum emission gives a flux den- 
sity of 1.56 mjy with a rms noise a = 0.75 mjy. We thus do 
not consider this a reliable determination of the flux density 
at 1 mm, and instead adopt a 3cr upper limit of 2.25 mjy. 

3.2 CO line emission 

Fig. [2] shows the channel maps for the CO(1-0) and CO(2- 
1) lines from the 60 km s _1 datacubes. To draw attention 
to the overall morphology and kinematics of the CO, we 
created integrated intensity maps, mean velocity field maps 
and position-velocity (PV) diagrams using the 40 km s _1 
datacubes, for better velocity resolution (Figs. [3] and |4]). 
To make these maps, we first created a smoothed cube by 



smoothing with a 2D gaussian (FWHM of 4 pixels) spatially 
and hanning smoothing by 3 channels in velocity. We then 
computed the moments by integrating in velocity over the 
pixels in the original cube that corresponded to pixels above 
3cr in the smoothed cube. As the emission is oriented primar- 
ily north-south, we created the PV diagrams by computing 
the zeroth moment of the masked cube in the Declination- 
velocity plane. This smoothing-masking procedure excludes 
much, but not all, noise from the moment maps. 

As shown in the integrated intensity maps of both tran- 
sitions, the CO is distributed somewhat asymmetrically, al- 
though with the signal just above 3cr in most of the velocity 
channels, we have to be wary of the contribution of noise 
even with the smoothing and masking procedure applied. 
The smaller CO(2-l) beam resolves two peaks in the CO in- 
tensity, one directly south of the galaxy centre and one to the 
north by north-east. The two clear CO(2-l) peaks suggest 
that the densest gas might be in a disc or ring seen nearly 
edge-on. The lower-resolution CO(1-0) data blend these two 
separate peaks into one peak elongated in the direction of 
the CO(2-l) peaks. In addition, the CO(1-0) data pick up an 
extended structure along PA~ 30°, especially prominent to 
the north-east. While this feature is absent in the CO (2-1) 
maps, it appears in both the 40 and 80 km s _1 channels of 
the CO(1-0) data, suggesting it is not just noise. The CO(2- 
1) observations must miss this molecular gas through some 
combination of the more limited primary beam size, noise, 
and a potentially lower excitation temperature of the in- 
falling gas. Overall, the molecular gas distribution suggests 
that we are observing the formation of an inner polar ring 
or disc, with gas in the outskirts of the inner kiloparsec still 
being accreted onto the settled structure. 

The velocity maps (Fig. [3} and PV diagrams (Fig. [4| 
support this interpretation. The overlay of the CO(1-0) PV 
contours on the CO(2-l) PV pixel map in Fig. |4]shows that 
the rotation of the main body of molecular gas is in good 
agreement for the two lines. Unfortunately, the quality of 
the data does not allow us to reliably determine whether 
the molecular gas is in a ring or a disc. For simplicity and 
given the hint of a flattening of the rotation at offsets larger 
than ~ 22", we refer to it as a disc for the remainder of 
the paper but refrain from making any claims that depend 
on whether it is a disc or a ring. The polar rotation is con- 
siste nt with the veloc ity field of the more extended ionized 
gas l|Sarzi et al. 2006). However, the CO(1-0) mean velocity 
map is irregular in the northern part of the galaxy. This ir- 
regularity is caused by the lower velocity (40 and 80 km s _1 ) 
molecular gas making up the extension to the north-east of 
the galaxy centre, as can be clearly seen in the CO (1-0) PV 
contours. The velocity and extent of this gas lead to the con- 
clusion that it must not yet be part of the central regularly 
rotating structure, but instead be in the process of accreting 
onto this structure. While other minor irregular features in 
the velocity maps may indicate further irregularities in the 
molecular gas kinematics, the accidental inclusion of a noise 
peak from a particular channel map could influence the ve- 
locity map in the same way, thus we hesitate to analyse these 
features. 

To construct integrated spectra and determine total 
fluxes, we selected the pixels in the central non-zero parts of 
the two integrated intensity maps as our spatial regions of 
interest. Integrating the original datacubes over these spa- 
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Figure 2. Channel maps of NGC 2768. Top: 12 CO(1-0) line. Bottom: 12 CO(2-l) line. The channels are 60 km s _1 wide and contours 
are plotted at -3 (dashed), 3, 6, and 9 <r with <r = 1.0 mjy beam -1 for CO(1-0) and a = 1.6 mjy beam -1 for CO(2-l). The number in 



the top left corner of each frame is the central velocity of that frame, relative to the observed central velocity of V sy 



1373 km s -1 . The 



clean ed beam is shown in the bottom left corner of each frame. The cross marks the centre of the 5 GHz radio continuum dNagar et al.l 
l2005fl . 



tial regions gave us the spectra sho wn in Fig. [5] Spectra f rom 
the 30m single-dish observations of lCombes et al.l |2007l ) are 
over-plotted in dotted lines. 

Comparison with the I Combes et al.l (120071 ) 30m data 
shows that the double-peaked profile seen in both lines in the 
interferometric data is absent in the 30m single-dish data. 
While the large amount of noise in the 30m spectra certainly 
contributes to the difference in the spectral profiles, we also 
note that both 30m spectra are biased toward the positive 
relative velocities seen on the northern side of the galaxy. If 
the 30m's pointing were off to the north, it would explain 
why intensity is preferentially missing at the negative veloc- 
ities. However, the intensities at positive relative velocities 
at 1 mm are also much lower than found with the interfer- 
ometer. This difference suggests that the pointing may have 
also been off in the east- west direction. The 1 mm spectrum 
suffered much more from this pointing error as the half- 



Table 2. CO fluxes and molecular mass estimates in NGC 2768. 



Line 


Instrument 


Flux 


M H2 






(Jy km s — 1 ) 


(10 7 M ) 


12 CO(1-0) 


30m 


10.4 ±1.6 ±1.0 


6.4 ± 1.0 ±0.6 


12 CO(1-0) 


PdB 


11.3 ±1.3 ±2.3 


6.8 ±0.8 ±1.4 


12 CO(2-l) 


30m 


18.6 ±2.5 ±1.9 




12 CO(2-l) 


PdB 


30.6±1.7±6.1 





Errors reported are random and systematic (flux calibration) 
errors, respectively. 



power beam- width is only 11". We note that errors in the 
flux calibrations could also explain some of the discrepancy. 

To estimate the total CO flux in each line, we inte- 
grated the spectra over the velocity range with observed 
emission, i.e. from 1193 to 1553 km s _ . Table [2] reports 
the calculated total fluxes, 11.3 Jy km s -1 at CO(1-0) and 
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Figure 3. Left: CO integrated intensity maps of NGC 2768. Right: CO mean velocity maps. Top: CO(l-O). Bottom: CO(2-l). The 
synthesized beam is shown in the bottom-left corner of each map. The black cross marks the centre of the 5 GHz radio continuum 
iNagar et al.ll2005l l . White pixels have been masked according to the procedure described in Section 3.2. 




Figure 4. Position- velocity diagram of NGC 2768 along the polar 
axis. The pixels show the CO(2-l) emission, while the contours 
show the CO(l-O) emission at 3, 6, and 9 xlO 8 Jy beam -1 km 
s . White pixels have been masked according to the procedure 
described in Section 3.2. 



30.6 Jy km s at CO(2-l) . Fluxes for the 30m observa- 
tions of lCombes et al. I (120071 ) are also reported; these fluxes 
have been calculated over identical velocity ranges as their 
interferometric counterparts. For CO (1-0), the fluxes agree 
within the noise and the errors in the flux calibration. For 
the CO(2-l) line, the 30m misses much of the flux seen by 
the interferometer, especially at negative relative velocities 
where the flux from the 30m is 8.2 Jy km s _1 less than for 
the interferometer. As we have noted above, we believe much 
of this difference to be due to a pointing error. 

We used the total fluxes in CO (1-0) to compute the to- 



tal molecular hydrogen mass using the formula M(H2) = 
(1.22 x 10 4 Mq) D 2 * Sco, where D is the distance mea- 
sured in Mpc and Sco is the total CO (1-0) flux. This for- 
mula comes from using the standard CO to H.2 conversion 
ratio 7V(H 2 )//(CO) = 3 x 10 20 cm' 2 , where N(U 2 ) is the 
column density of Eh and I(CO) is the CO(1-0) intensity in 
Kkms" . The total H 2 mass detected by our interferometric 
observations is then 6.8 x 10 7 Mq. As nearly all the CO(1-0) 
flux comes from the nascent polar disc seen in the integrated 
intensity maps, this gives a rough estimate of the molecular 
mass already in the central polar disc. We estimate the in- 
tegrated line ratio R 2 i = CO(2-1)/CO(1-0)^ 0.7, indicating 
optically thick molecular gas that is subthermally excited, as 
expected for a non-starburst galaxy. This ratio is computed 
assuming we have mapped the full extent of the CO in both 
lines in our interferometric maps. As the CO distribution in 
both lines is only about half the diameter of the HPBW of 
the primary beam, this assumption is reasonable. 



4 DISCUSSION 

4.1 Origin of the Molecular Gas 

With H I data from iMorgariti et al l i|2006h . dust maps from 
HST images (e.g. Fig. [6}, and the molecular gas maps of 
section 3, we now have a fairly complete census of the cool 
ISM in and around NGC 2768. The most striking feature, 
besides the fact that the ISM lies in the polar plane, is the 



6 A. F. Crocker et al. 



0.05 - 




-500 500 

Relative Velocity (km s ) 




-200 200 

Relative Velocity (km s ') 

Figure 5. Solid lines show integrated spectra of the CO(1-0) line 
(top) and CO(2-l) line (botto m) in NGC 2768. The dotted lines 
show the single-dish spectra of lCombes et al" | ^0071) . 



clear link between all components. All observed phases of 
the cool ISM extend to the north-east, first in CO, then 
dust, then HI (Fig. [B). This shared asymmetry ranges from 
the inner kiloparsec in the molecular gas to over 30 kpc in 
H I. This link between the different ISM components clearly 
indicates that they were all accreted from the same source. 
It also suggests that the gas accretion will continue for some 
time, perhaps in discrete episodes. 

This extended and asymmetric morphology along with 
the misaligned kinematics of the cool gas also firmly indi- 
cates an external origin. We consider three options: a minor 
merger with a gas-rich dwarf, tidal accretion from a gas-rich 
galaxy, and accretion from an HI filament. The regularity of 
both the photometry and stellar kinematics of NGC 2768 ar- 
gues against the recent merger of a gas-rich dwarf galaxy. No 
stellar tidal tail nor other irregularity has been observed, al- 
though the data we have access to do not pose particularly 
stringent upper limits (i.e. the optical photometry is rela- 
tively shallow). The accretion of primordial (metal- free) HI 
scenario fails to explain both the dust and CO observed. We 
thus consider tidal accretion as the most likely scenario. 

The loose-group centred on N GC 2768, cata logued as 
Lyon Group of Galaxies (LGG) 167 (|Garcialll993h . contains 



many potential gas-d onor galaxies (see Ta ble [3). The H I 
datacube obtained bv lMorganti et al] (|2006l ) reveals neutral 
hydrogen around 4 lower- luminosity galaxies within a 150 
kpc projected separation from NGC 2768. The H I veloc- 
ities indicate group membership for all of these galaxies. 
Considering the wide range (0.1 to 6) of H I mass to blue 
luminosit y ratio (Mu i / Lb) values found in late- type dwarf 
galaxies (|Swatersl 1 1999T ) , any of these galaxies could have 
provided the ~ 10 s M of H I found around NGC 2768 
without requiring an unreasonable amount of H I for its blue 
luminosity. 

Including all galaxies within 400 kpc in projected sepa- 
ration and 400 km s _1 in relative line-of-sight velocity from 
NGC 2768, we created a list of known galaxies that may 
have recently interacted with NGC 2768 (Table IS}. Given 
the r elaxed appearanc e of the ionised gas in the SAURON 
map IjSarzi et al.ll2006f ). we estimate that this gas must have 
been present for at least 3-10 dynamical timescales, or 0.2- 
0.7 Gyr. Assuming a gTOUp = 200 km s _1 as a reasonable 
three-dimensional velocity dispersion for a small group, we 
estimated the time since each galaxy could have interacted 
with NGC 2768, 



Y w group i-ivsys 

where AR is the projected separation and AV sys the line-of- 
sight velocity relative to NGC 2768. Values of this estimated 
interaction time are listed in Table [3] This analysis shows 
that the nearby Scd galaxy UGC 4808 is most likely to have 
interacted with NGC 2768 recently and provided the cool 
ISM. UGC 4808 can be seen in the bottom right corner of 
Fig. [6] Its alignment with the large-scale H I supports this 
conclusion. 

4.2 Fate of the Molecular Gas: Star Formation? 

In NGC 2768, we observe no definitive signs of young stars 
as seen in the other five SAURON early-type galaxies with 
mapped CO emission (NGC 2865, Schinnerer & Scoville 
2002; NGC 3032, NGC 4150 NG C 4459 and NGC 4526, 
lYoung. Bureau fc Cappeliaril 120081 ). Excluding NGC 2685, 
in which the CO lies outside of the SAURON field-of- 
view, young stars in these galaxies are traced by high H/3 
linestrengths (all), cold stellar components that share the 
CO kinematics (all except NGC 3032), and low [OKI] to H/3 
emission line ratios (all except NGC 4150). 

The H/3 linestrength index of NGC 2768 is 1.70 A 
within a central 9" radius aperture (one eighth of the ef- 
fective radius) and remains near this level over the entire 
SAURON field (Kuntschner et al. 2006; Fig. |7J. Higher 
resolution but lower signal-to-noise data from the OASIS 
integral-field spe ctrograph give slightly larger values for the 
H/3 linestrength dMcDermid et all 2006). However, the dis- 
crepancy between the OASIS and SAURON values is prob- 
ably attributable to the difficulty of subtracting the strong 
H/3 emission before fitting the linestrengths, and the higher 
signal-to-noise data from SAURON provide a better esti- 
mate. Combining line-strength indices (H/3, Fe5015, and Mg 
b) with stellar population models, Kuntschner et al. (in 
preparation) show that NGC 2768 has a fairly uniform, old 
stellar population. However, linestrengths are not sensitive 
to very recent star formation nor to a very small percentage 
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Table 3. Basic properties of galaxies near NGC 2768. This table lists all galaxies with a projected distance differing by less than 400 kpc 
and a redshift differing by less than 400 km s — 1 from NGC 2768. Only galaxies with known redshifts are included. 



Name 


A8 


AR 


AV SyS 


B T 


Mhi 


M H i /L b 


At 


Type 






(kpc) 


(km s- 1 ) 


(mag) 


(10 8 M ) 


(M /L ) 


(Gyr) 


(NED) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 2768 








10.60 


1.7 


0.007 




E 


UGC 04808 


11.'5 


72 


-78 


15.03 


1.6 


0.41 


0.38 


Scd 


PGC 025982 


16.'5 


103 


-236 


14.71 










MAILYAN 038 


19.'8 


121 


187 




0.3 




1.71 




PGC 2599651 


21'.8 


136 


197 


17.05 


0.4 


0.63 


3.83 




PGC 2601563 


22.'5 


111 


-325 


17.36 


0.6 


1.4 






SDSS J091531.98+594948.5 


31f9 


200 


-5 


17.14 






0.98 




NGC 2742 


40.'2 


252 


-84 


11.38 


18.4 


0.16 


1.36 


SA(s)c 


NGC 2726 


50.'6 


317 


115 


12.56 


1.8 


0.05 


2.25 


Sa? 



(1) Galaxy name. (2) Angular separation. (3) Projected separation (assuming a distance of 21.5 Mpc). (4) Relative line-of-sight 
velocity. (5) Total corrected apparent blue magnitude (HyperLEDA). (6) Mass of neutral hydrogen (Morganti et al. 2006, except for 
NGC 2742 - Broeils & van Woerden 1994 and NGC 2726 - Haynes et al. 1988). (7) HI mass to total blue luminosity ratio. (8) 
Estimate for time since a possible interaction with NGC 2768. (9) Galaxy classification (NED). 
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Figure 6. The cool ISM in NGC 2768. The left panel shows the H I surface density distribution in contours of 2, 4, 8, 16 and 
32 X 10 19 atom cm 2 plotted over an R-band Digital Sky Survey (DSS) image. UGC 4808 is visible in the bottom right-hand corner. 
The right panel shows the CO(1-0) intensity distribution in contours of 0.5, 0.75, 1.0 and 1.25 Jy km s -1 plotted over a dust map of 
NGC 2768 obtained by dividing the HST F814W ACS image by an ellipse model fit. In the right panel, the outer circle indicates the 
44" HPBW of the CO(1-0) primary beam of the interferometer and the arrows indicate the galaxy centre. The black ellipses in the 
bottom left corners of both panels show the synthesized beams of the H I and CO observations, respectively. 



of young stars superposed on a dominant older stellar pop- 
ulation, so linestrengths cannot rule out particularly recent 
and/or weak star formation. 

The SAURON map of the stellar velocity dispersion 
shows a dramatic drop of ~ 45 km s~ in the innermost 
8" (Emsellem et al. 2004; Fig. [7J). This drop in dispersion is 
also clearly seen in lo ng-slit data l|Simien fc Prugniel|[l997l ) 
and in OASIS data (|McDermid et al.1 120061 ). Young stars 
born from dynamically cold gas in a cicumnuclear disc can 
produce such a dro p in dispersion ijEmsellem et al.l l200ll : 
IWozniak et"al1l2003h . but they should also produce a notice- 
able effect in the velocity map. However, the velocity field of 



NGC 2768 lacks any signature of a dynamically cold stellar 
disc or ring. Thus the stellar velocity dispersion drop cannot 
be directly attributed to the formation of a stellar compo- 
nent either from the molecular polar ring itself, or from gas 
settled into the main pla ne of the galaxy as may occu r during 
polar accretion events (|Bournaud fc C ombes 2003 ; ). Addi- 
tionally, it is difficult to associate this drop to a young stellar 
population, as no coincident increase in H/3 linestrength is 
observed. We conclude that no kinematic component visible 
in the SAURON maps is directly related to the molecular 
gas. 

Particularly low ratios of [O III] to H/3 [log([0 III] /H/3) 
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Figure 7. Top Left: CO(l-O) contours of NGC 2768 over the SAURON H/3 absorption linestrength map l lKuntschner et al-lliou^) . No 
signif icant change in lincstrength is coincident with the molecular gas. Top Right: CO(2-l) contours over the stellar velocity dispersion 
map l lEmsellem et alj|2004). The CO is co-spatial with the dispersion drop. Bottom Left: CO(l-O) contours over a map of the equivalent 
width of H/3 emission l|Sarzi et al1l2006h . Note the increase in equivalent widt h near the extend ed CO peak to the north-east. Bottom 
Right: CO(l-O) contours over a map of the flux ratio of [O III] to H/3 emission llSarzi et alj|200rj) . 



< —0 .2] indicate st ar formation as the source of ionisa- 
tion (|Ho et al.lll997h . In NGC 2768, log([0 III]/H/3) ranges 
from about 0.0 to 0.5. In this range, the ionisation can 
come from sources other than star formation. The four 
additional possibilities are heating by an X-ray halo, an 
AGN shocks or post-AGB stars. NGC 2768 lacks an X-ray 
halo l|Komossa et al.lll999l ). so this option is ruled out. The 
H/3 flux in NGC 2768 does not drop off with the square 
of the radius as it should if io nised solely by an AGN 
dHatch, Crawford fc Fabianll2007l ). therefore everything out- 
side the central 2" must be ionised by another source (Sarzi 
et al., in preparation). 

The most prominent feature of the ionised gas in 
NGC 2768 is a roughly spiral-sha ped region of inc reased 
[O III] equivalent width (see fig. 4b. lSarzi et~ai1l2006l ). This 
strong enhancement in [OIII] without a corresponding strong 
enhancement in H/3 is seen in other SAURON galaxies 
(NGC 2974, NGC 3414) and must be caused by shock ion- 
isation of the gas, as post-AGB stars are too symmetric to 
cause such asymmetric ionisation and star formation would 
not cause such an increase in [OIII]. Yet the smooth and cir- 
cular H/3 emission suggests that an ionisation source other 
than shocks is responsible for the ionisation outside of the 
spiral pattern. Both star formation and post-AGB stars re- 
main options. Post-AGB stars must be present in such an old 
galaxy and thus must play a role. However, the H/3 equiva- 
lent width map from SAURON shows a weak enhancement 
along the extended molecular gas seen to the northeast of 
the galaxy centre (see Fig. [7}, pointing to star formation in 
this ionised gas. 



The FIR luminosity (Lfir) to molecular mass (Mh 2 ) 
ratio also provides possible evidence of star formation 
jCombes et al.ll2007l ). We computed a bolometric FIR flux 
for NGC 2768, 

FIR = 1.26 x 10 -14 Wm- 2 (2.58 S 6 o M m + SWm), 

where Ssoum and Sioo/jm are IRAS fluxes in Jy. The total 
FIR luminosity is then L FIR = (FIR) AttD 2 = 4.3 x 10 8 
Lq , giving a Lfir to Mh 2 ratio of 6.3 L(n/M(Ti, a reasonabl e 
value for a star- forming spiral l|Devereux fc Hameed|[l997l ). 
In spiral galaxies, FIR emission traces star formation activ- 
ity and the ratio of Lfir to Mh 2 is assumed to measure how 
efficiently stars are being formed out of the molecular hydro- 
gen. If star formation is the dominant process contributing 
to FIR emission in NGC 2768, then star formation is pro- 
ceeding at a normal efficiency. However, other sources (such 
as an AGN) may significantly contribute to the FIR emis- 
sion in early-type galaxies, weakening this conclusion. We 
can still use the FIR luminosity to give an upper limit to 
the star formatio n rate, deriving a n upper limit of 0.07 solar 
masses per year (|Kennicuttlll9981 . eqn. 4). Thus if star for- 
mation is present in NGC 2768, it is occurring at a very low 
rate, and the timescale to exhaust the supply of molecular 
gas is around 1 Gyr. 

As we have only unclear evidence of ongoing star for- 
mation in NGC 2768, we consider whether the molecular 
gas disc might be gravitationally stab le and thus unable to 
form stars. Theoretical models from iKawata et al.l (|2007t ) 
suggest that circumnuclear discs in massive galaxies are 
more stable than those in less massive galaxies. A stabil- 
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ity analysis of the molecular circu mnuclear disc in 3C 31 
based on Toomre's Q-parameter (jToomre! 1 1964 ) suggests 
that i t is gravitationally stable and thus unable to form 
stars (|Okuda et al.| [2005). Our data are insufficient to cal- 
culate Toomre's Q parameter by thems elves, but comparing 
with the models of lKawata et al.l l|2007h (see their fig. 5) , we 
find that the 0.5 kpc radius and 6.8 x 10 7 Mq mass of the 
molecular disc in NGC 2768 should be unstable, with a Q 
parameter less than 1. Thus it is reasonable to expect star 
formation in the polar disc of NGC 2768. 



4.3 Origin of the mm continuum 

From radio and X-ray observations, we know that NGC 2768 
is a low- luminosity AGN. Very Large Array (VLA) and 
Very Long Baseline Array (VLBA) observatio ns indicate a 
comp act radio source with no evidence of jets l|Nagar et al.l 
2005). Most if not all X -ray emission is consistent with a 
point-source origin, and the X-ray count rate differed be- 
tween observations a month apart, indic ating that the X-ray 
source is variable (jKomossa et alj fl999). 

The few LLAGN that have been studied at millimetre 
frequencies show a range of spectral shapes attributed to 
varying contributions from optically thin synchrotron emis- 
sion, self-a bsorbed synchrotron emission, and free-free emis- 
sion (see INagar et al] |2002| ; iDoi et alj 120051 ; iKrips et all 
120071 ). In Fig. |8] we plot the radio through FIR SED 
of NGC 2768, which is similar to those of some of the 
other LLAGN previously observed. Values and references are 
listed in Table [4] We note that the flux densities at different 
frequencies are necessarily at different resolutions and thus 
are not entirely comparable. However, as no extended struc- 
tures are observed at the milliarcsecond resolution of the 
VLBA at 5 GHz, and the 1.4, 15, and 115 GHz detections 
are all consistent with a point-source origin, we assume all 
these flux densities are dominated by an unresolved source 
and can thus be compared without correction. The upper- 
limit at 230 GHz is also for a point source at the phase centre 
at the resolution of 1'.'2 x 1'.'2, while the 240 GHz detection is 
a bolometric measurement from an 11" beam, which could 
inclu de a contribution from cold dust (|Wiklind fc Henkel 
1 19951 ) as well as from a central AGN. These two observa- 
tions are thus not as contradictory as they appear in the 
SED plot. However, we note that our map of 1 mm contin- 
uum emission does not trace the dust pattern seen in Fig. [5] 
nor do we see evidence of more emission in the central 11" af- 
ter smoothing our data to this resolution. 

First, the lack of strong, coincident 1 mm continuum 
emission rules out cold dust as the origin of the 3 mm con- 
tinuum emission, as the thermal spectrum from dust drops 
sharply towards mm wavelengths (see the dust curves in 
Fig. [SJ. We can also rule out free- free emission using the 
formula in IDoi et al.1 (|2005T ) which assumes an electron tem- 
perature of 10 4 K and N(He+)/N(H+)= 0.08: 



V mJ y/ 



3.57 x 10 1 



F(H0) 



erg cm - 



(ghz) 



We use the tota l flux in H/3 from the SAURON data 
l|Sarzi et al.ll2006l ). which gives an overestimate for the cen- 
tral radio emission, as most of the H/3 emission comes from 
ionising sources other than star formation and only the flux 



Table 4. Spectral energy distribution of NGC 2768. 



Instrument Frequency (GHz) Flux Density (mjy) Rcf. 



WRST 


1.4 


10.9 


1 


VLBA 


5 


7.3 


2 


VLA 


15 


8.2 


2 


PdB 


115 


1.94 


3 


PdB 


230 


< 2.25 


3 


IRAM 30m 


240 


10.1 


1 


MIPS 


1875 


377 


5 


IRAS 


3000 


1370 


6 


MIPS 
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5 


IRAS 
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Figure 8. Spectral energy distribution of NGC 2768 from the 
radio to the far infrared (FIR). The solid line is an a = —0.5 
power law fit, modelling synchrotron emission. The dot-dashed 
line estimates the maximum contribution from free-free emission. 
The dashed curve extrapolates from the FIR dust bump to mm 
wavelengths using a f3 = 1 emissivity law; the dotted curves does 
the same with /3 = 2. 



in the central pixels should correspond to H II regions that 
could appear as a central radio point source. Nevertheless, 
Fig. H] reveals that, even with this overestimate, free-free 
emission is far too faint to account for the 3 mm continuum. 

Having ruled out dust and free-free emission, processes 
directly related with an accreting supermassive black hole 
are likely to be the source of the 3 mm continuum detected in 
NGC 2768. Optically thin synchrotron emission from com- 
pact jets or other ejecta will produce a spectrum with a 
spectral index of a ^ —0.5. The logarithmic least-squares 
fit of an a — —0.5 power law is shown in Fig. [8] It roughly 
fits the 3 radio data points and the 3 mm continuum emis- 
sion, but the low value at 1.4 GHz and the high value at 
15 GHz hint at a flatter or inverted spectrum, suggesting 
synchrotron self-absorption plays a role. 
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5 CONCLUSIONS 

The interferometric observations of the CO (1-0) and CO (2- 
1) lines presented here reveal a kiloparsec polar ring or disc 
of molecular gas within the early-type galaxy NGC 2768. 
This molecular gas reveals the potential of this galaxy 
to form kinematic substructure through star formation. 
However, neither the kinematics nor the linestrengths of 
NGC 2768 as measured by SAURON show clear evidence 
of a young stellar population. This limits any star formation 
coincident with the molecular gas to be very recent and/or 
too small a luminosity-weighted fraction of the background 
old stellar population to contribute to the linestrengths 
and kinematics. The asymmetric and unrelaxed nature of 
the CO (1-0) emission suggests that the molecular gas in 
NGC 2768 was recently accreted and is possibly still be- 
ing accreted, in agreement with the much more extended 
but disturbed HI distribution. Searching the neighbourhood 
of NGC 2768 for potential gas-donor galaxies, many H I- 
rich low-luminosity galaxies are found, all of which could be 
the source of H I seen in and around NGC 2768, although 
UGC 4808 seems the most likely, being closest in projected 
distance and aligned with the extended H I distribution. 

These CO observations of NGC 2768 probe new regions 
of parameter space in terms of molecular gas in early-type 
galaxies, with recently accreted, polar-rotating gas of only 
6.8 x 10 7 M . As surveys have indicated lower CO detec- 
tion rates in higher-luminosity E an d SO galaxies l|Lees et al.l 
Il99ll ; ISage, Welch fc Yound [20071 ) . NGC 2768 is also no- 
table fo r being the brightes t in B-band of the CO detec- 
tions in lCombes et al.l (|2007T l. Whether the lack of observed 
star formation has to do with the (short) time since the gas 
has been accreted, the (low) mass of molecular gas avail- 
able for star formation, the polar dynamics of the gas, or 
some host-galaxy property will take more investigation to 
discover, including an analysis of the molecular gas distri- 
butions of other early-type galaxies with integral-field data. 
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